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ABSTRACT 

Comparing models of Simple Stellar Populations (SSP) with observed line strengths generally provides a tool 
to break the age-metallicity degeneracy in elliptical galaxies. Due to the wide range of Balmer line strengths 
observed, ellipticals have been interpreted to exhibit an appreciable scatter in age. In this paper, we analyze 
Composite Stellar Population models with a simple mix of an old metal-rich and an old metal-poor component. 
We show that these models simultaneously produce strong Balmer lines and strong metallic lines without invoking 
a young population. The key to this result is that our models are based on SSPs that better match the steep increase 
of H/3 in metal-poor globular clusters than models in the literature. Hence, the scatter of H/3 observed in cluster 
and luminous field elliptical galaxies can be explained by a spread in the metallicity of old stellar populations. We 
check our model with respect to the so-called G-dwarf problem in ellipticals. For a galaxy subsample covering 
a large range in ISOO-V colors we demonstrate that the addition of an old metal-poor subcomponent does not 
invalidate other observational constraints like colors and the flux in the mid-UV. 

Subject headings: galaxies: elliptical and lenticular, cD - galaxies: stellar content - galaxies: formation - 
galaxies: abundances - galaxies: fundamental parameters (ages, metallicities) 



1. INTRODUCTION 

More than 20 years ago it has been recognized that the mod- 
eling of the spectral energy distribution of ellipticals is affected 
by an ambiguity in age and metallicity (Faber 1972; O'Connell 
1976), which has turned out to be a general complication in 
population synthesis (Renzini 1986). However, considering 
Simple Stellar Populations (SSP) in the two-parameter space of 
Balmer and metallic lines, the age-metallicity degeneracy can 
be broken (Gonzalez 1993; Worthey 1994). The major reason 
for the success of this strategy is that the Balmer line strengths 
of SSPs are predominantly age sensitive at metallicities above 
^ 1/3 Zq that are supposed to be the only relevant for ellip- 
tical galaxies. Strong H/3 lines are thus taken as an indication 
for young (intermediate-age) populations, the observed scatter 
is interpreted as a considerable spread in age (e.g., Faber et al. 
1995). 

An alternative path to obtain blue stars and hence strong 
Balmer lines is to consider old metaZ-/?oor populations. In this 
paper we follow this approach and compute composite stellar 
populations that contain an old metal-poor subcomponent. The 
principal focus is to check if a combination of only old pop- 
ulations can reproduce strong H/3 without invalidating further 
constraints for ellipticals like metallic indices, colors and spec- 
tral energy distributions. 

The paper is organized as follows. In Section ^ we cali- 
brate the SSP model indices and spectral energy distributions 
on globular cluster data. The composite models and their appli- 
cation to galaxy data are presented in Section ^ In Sections Q 
and ^ we discuss and summarize the results. 

2. SIMPLE STELLAR POPULATION MODELS 

We compute a new set of SSP models based on the popula- 
tion synthesis presented in Maraston (1998), in which the fuel 
consumption theorem (Renzini & Buzzoni 1986) is adopted to 
evaluate the energetics of the post main sequence phases. The 
new SSP models cover the metallicities -2.25 < [Fe/H] < 0.5 



and ages 10^-2- 10'" yr The input stellar tracks are taken 
from Bono et al. (1997) and S. Cassisi (1999, private commu- 
nication). These new models will be discussed in detail in a 
future paper (C. Maraston, in preparation). The calibration of 
the model colors on globular clusters is presented in Maraston 
(1998, 2000). Here we show the calibration of indices and spec- 
tral energy distributions of the SSP models that are relevant for 
the present study. 



2.1. Spectral line indices 

Synthetic line indices for SSPs are computed using the fit- 
ting functions from Worthey et al. (1994). Fig. |l] shows H/3, 
Fe5335, and Mgb as a function of [Fe/H] for galactic globular 
clusters from various data sets (see the caption). Our SSP mod- 
els of fixed age (15 Gyr) and various metallicities are plotted as 
solid lines. Dotted lines are the models from Worthey (1994) 
forf = 17Gyr. 

As metal-poor stars ([Fe/H] ^ -1.3) spend their horizontal 
branch lifetimes at temperatures Tgff ;> 5500 K, H/3 is steeply 
rising with decreasing [Fe/H] (upper panel in Fig. nl). This ob- 
served feature is the key for the present analysis and is well 
reproduced by our SSPs. The models from Worthey (1994), in- 
stead, predict a shallower trend, which implies that the very 
metal-poor globulars are younger than the more metal-rich 
ones. The synthetic Fe5335 indices are in good agreement 
with the Covino, Galletti & Pasinetti (1995) data, while for 
[Fe/H] > -1.5 Burstein et al. (1984) and Trager (1998) derive 
lower values. Mgb indices seem slightly overestimated by the 
models at low metallicities. 

Our synthetic metallic indices Fe5335, Mgb and the Balmer 
line H/3 for [Fe/H] ^ -1.5 agree well with the models from 
Worthey (1994), because both sets of models are based on the 
same fitting functions (Worthey et al. 1994). The uncertainty 
introduced by the use of different fitting functions is discussed 
in Maraston, Greggio & Thomas (1999) and L. Greggio & 
C. Maraston (in preparation). The low H/3 values given by 
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[Fe, H] 

Fig. 1. — Calibration of the SSP models on galactic globular clusters. Data are from Burstein et al. 1984 (squares), Covino, Galletti 
& Pasinetti 1995 (circles), and Trager 1998 (triangles). [Fe/H] (Zinn & West (1984) scale) is taken from Harris (1996). Solid lines 
are our SSPs for age f = 15 Gyr. Worthey (1994) SSPs for t = 11 Gyr are shown as dotted lines. 



Worthey (1994) at very low metallicities are likely due to the 
assumption of too red horizontal branch morphologies in the 
models. 

As discussed in Worthey et al. (1994), iron indices are fa- 
vorable as metallicity indicators, because the analysis with Mg- 
indices is severely affected by a-enhancement (Tripicco & Bell 
1995): SSP models based on solar abundance ratios predict - 
for a given Fe index - weaker Mg indices than measured in 
elliptical galaxies (Worthey, Faber & Gonzalez 1992; Greggio 
1997). As at high metallicities the fitting functions of Fe5270 
are significantly more uncertain (Maraston et al. 1999), in the 
following analysis we focus on H/3 and Fe5335. We addition- 
ally consider the color B-K as metallicity indicator which turns 
out to be in good agreement with Fe5335, while Mgb system- 
atically implies higher metallicities (see Section ||). 

Combining H/3 with Fe5335, for 47 Tuc we obtain a 'spec- 
troscopic age' of 15 Gyr (Fig. 0) in agreement with the age 
derived from the color-magnitude diagram (14 ± 1 Gyr, Richer 
et al. 1996). This example shows that the H/3-Fe5335 plane 
is successful in breaking the age-metallicity degeneracy, when 



applied to simple systems like globular clusters. Note that the 
high resolution Balmer index H7HR, instead, leads to a 'spec- 
troscopic age' in excess of 20 Gyr (Gibson et al. 2000). 

2.2. Spectral energy distributions 

The evolutionary synthesis code of Maraston (1998) is up- 
dated for the computation of the SSP spectral energy distri- 
butions (SEDs). The spectral library of Lejeune, Cuisinier & 
Buser (1998) is adopted to describe the stellar spectra as func- 
tions of gravity, temperature and metallicity. 

The tightest constraint on the amount and the metallicity of 
metal-poor stars in elliptical galaxies comes from the flux in 
the mid-UV (2000-4000 A). In Fig. | we present the calibra- 
tion of synthetic metal-poor SEDs on galactic globular cluster 
data in the wavelength range A = 1500-3300 A (van Albada, 
de Boer & Dickens 1981). The model spectra are in excellent 
agreement with the data in the whole metallicity range. The in- 
crease of the flux shortward 3000 A with decreasing metallicity 
is due to a hotter main sequence turn-off and a bluer horizontal 
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Fig. 2. — Calibration of the synthetic SEDs in the UV on galactic globular clusters. Data (squares) are taken from van Albada, de 
Boer & Dickens (1981), [Fe/H] (Zinn & West (1984) scale) is adopted from Harris (1996). The filled squares (central panel) show 
Ml 3 that has an extreme blue horizontal branch morphology. The solid lines are our model SEDs (f = 15 Gyr) with metallicities as 
indicated in the panels. 



branch. Note, however, that there are globulars of intermediate 
metallicity ([Fe/H]~ -1 .5) that show horizontal branches bluer 
than what it is expected with the canonical mass-loss on the Red 
Giant Branch (van Albada et al. 1981). This case is shown in 
the central panel of Fig. |for the globulars M13 (NGC 6205), 
filled squares) and NGC 5272, which have the same metallic- 
ity [Fe/H] -1.54 but different horizontal branch morpholo- 
gies. M13 contains extreme blue horizontal branch stars ('blue 
tails'), and ex hibits an excess of flux shortward ~ 2000 A (see 
Section |3.2.l|). 



As blue tails are preferentially found in denser and more con- 
centrated globulars, it is likely that the particular dynamical 
conditions in such 'second parameter' globular clusters lead to 
enhanced mass-loss and extremely hot horizontal branch stars 
(Fusi Pecci et al. 1993), which may not apply to the average 
stellar field population in galaxies. Since in early-type galaxies 
globular clusters contribute less than 0.6 per cent to the total V- 
light (Ashman & Zepf 1998), we adopt the canonical mass loss 
ry = 0.33 for [Fe/H]~ -1 .5 calibrated on 'normal' globulars like 
NGC 5272 (Fusi Pecci & Renzini 1975). 

3. RESULTS 

In the following we compare old (> 12 Gyr) composite 
stellar populations with data of elliptical galaxies. The mod- 
els consist of two components: a major metal-rich population 
(Z = Znjaj) and a metal-poor subpopulation (Z = Zjub)- More 
complex models with a continuum of metallicities will be sub- 
ject of a future paper. 



3.1. Indices and colors 

Fig. H displays H/3 vs. Fe5335 (left-hand panel) and H/3 vs. 
B-K (right-hand panel) for composite models in which the ages 
of the major population and of the metal-poor subcomponent 
are fixed to 14 Gyr and 15 Gyr, respectively. The contribution 
from the low-metallicity population is 10 per cent by mass. 

The grid in Fig. ^ shows models for various metallicities of 
the two components. Solid fines connect models with constant 



metallicity of the major component Z^aj = 0.5, 1,2,3 Zq (from 
left to right). The metallicity of the subcomponent is Zsub = 0.05 
(large open squares) and Zsub = 0.005 Zq (large open circles). 
For comparison, simple stellar populations with metallicities 
Z = Zmaj and t = f^aj are shown as large open triangles. The 
dotted (horizontal) lines mark models of constant Zsub- 

The perturbation of metal-rich simple stellar populations 
with metal-poor stars leads to lower Fe5335 line strengths, 
bluer B-K colors, and stronger Balmer lines. The models 
with Zsub = 0.005 Zq reach H/3 2 A. Observed line strengths 
(Gonzalez 1993; Kuntschner & Davies 1998; Mehlert et al. 
2000) and colors (Pahre 1998) of elliptical galaxies and lentic- 
ular galaxies are shown as small filled and small open symbols, 
respectively. The present set of composite models match the 
area that is covered by the majority of the galaxy data. Note that 
the objects exhibiting H/3 > 2 A are either classified as lenticu- 
lar galaxies (small open symbols) or they are intermediate-mass 
^eWeflipticals (NGC 1700, NGC 3377, NGC 5831, NGC 6702, 
NGC 7454) or dwarf eflipticals (M32 (NGC 221), NGC 4489) 
from the Gonzalez (1993) sample. All luminous field ellipti- 
cals and all cluster ellipticals (Coma, Fornax, and Virgo) in 
Fig. ^ instead, can be modeled with the old composite popula- 
tions introduced above. It should be emphasized that the bulk 
of ellipticals scatter around H/3 ^ 1.6 A, so that less than 10 
per cent contribution from metal-poor stars is required in most 
cases. The observed range H/3 ~ 1.3-2.0 is then due to a spread 
in the metallicity and/or the weight of the metal-poor subcom- 
ponent. The analysis with simple stellar population models, 
instead, yields an a^e-spread of 7 - 15 Gyr. 

Another striking feature of Fig. ^ is that the positions of the 
galaxy data relative to each other are the same in the Fe5335 and 
the B-K diagrams. Moreover, the locations of the model grids 
relative to the data points in the two panels are in good agree- 
ment, so that both the metallic index Fe5335 and the color B-ZT 
consistently constrain the metallicity range for the major com- 
ponent to 1 - 3 Z0. This strongly reinforces the use of Fe5335 
as metallicity indicator (see Section H). In Fig. Qwe show the 
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Fig. 3. — H/3 vs. Fe5335 and B-K for elliptical (small filled symbols) and lenticular (small open symbols) galaxies. Data are from 
Gonzalez 1993 (circles, Virgo and field, /?e/8 aperture, labels indicate NGC numbers), Kuntschner & Davies 1998 (squares, Fornax), 
and Mehlert et al. 2000 (triangles. Coma). B-K colors are taken from Pahre (1998). The error bars denote average Icr errors. SSP 
models for Zssp = 0.5, 1,2,3 are large open triangles, with Zssp increasing from left to right. Composite models for Z^^j = Zssp 
and Zsub = 0.05,0.005 Zq are indicated by large open squares and large open circles. Solid and dotted lines denote fixed Zmaj and 
fixed Zsub, respectively. The metal-poor subpopulation contributes 10 per cent by mass. The ages of the major component and the 
subcomponent ai^e 14 Gyr and 15 Gyr, respectively. 



same data and the same model grid in the H/3-Mgb plane. As 
discussed in Section 2.1, Mg lines of ellipticals are generally 
stronger than any SSP model. This feature is attributed to an 
enhancement of Mg in luminous ellipticals. As a consequence 
the metallicities derived from the Mg-indices are substantially 
higher and provoke severe inconsistencies with observed colors 
(see also Saglia et al. 2000). As the latter are much less affected 
by abundance ratio anomalies, the HjS-B-K plane represents the 
more reliable constraint on stellar population models. 

Note that the H/3 line index additionally suffers from un- 
certainties in the correction for contamination by emission 
(Gonzalez 1993), which can lead to an underestimation of 
the H/3 line strength. This effect may explain the low values 
H^ ^ 1 .5 A that fall below the model grid in Fig. |[ 

3.2. Constraints from the mid-UV 

The fluxes radiated by metal-poor stellar populations are high 
in the mid-UV (2000-4000 A), and decrease shortwai'd 2000 A 
(Fig. The so-called UV-upturn below 2000 A observed in 
elliptical galaxies (Burstein et al. 1988), instead, is most likely 
due to late evolutionary phases of old metal-rich populations 
(Greggio & Renzini 1990; Dorman, O'Connell & Rood 1995; 
Yi, Demarque & Oemler 1998; Greggio & Renzini 1999). The 
amount of metal-poor stars in elliptical galaxies is therefore 
tightly constrained by the flux in the mid-UV around 2500 A. 



A model based on the closed-box metallicity distribution is 
not compatible with the mid-UV flux of the elliptical galaxies 
NGC 4649 and NGC 1404 (Bressan, Chiosi & Fagotto 1994) 
and of the M31 bulge (Worthey, Dorman & Jones 1996), which 
is referred to as the G-dwarf problem for ellipticals. On the 
other hand, Lotz, Ferguson & Bohlin (2000) find that the spec- 
tra of ellipticals in the mid-UV even require the addition of a 
metal-poor component. In the following we investigate if the 
proportions of metal-poor stars required to explain the H/3 line 
strengths are consistent with the fluxes observed in the mid-UV. 

The objects with the strongest H/3 (^ 1.65) that we ob- 
tain from cross correlating the samples of Gonzalez (1993), 
Kuntschner & Davies (1998), and Mehlert et al. (2000) with 
the Burstein et al. (1988) sample ('quiescent objects') are M31 
(NGC 224), NGC 4472, and NGC 3379. We additionally con- 
sider NGC 4649, which exhibits a low W(3 = 1 .40. These objects 
cover a large range in 1500-y colors, NGC 4649 having the 
strongest UV-upturn. In the following we check if the compos- 
ite models for these specific galaxies are in agreement with their 
SEDs. In order to constrain the flux from the metal-poor pop- 
ulation in the mid-UV, it is necessary to model the UV-upturn. 
Following the description by Greggio & Renzini (1990), we in- 
clude hot (r 25 , 000 K) late evolutionary phases in the metal- 
rich components, such that the observed UV rising branches are 
reproduced. 
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Fig. 4. — H/? vs. Mgb. Symbols and linestyles like in Fig. || 



The main model parameters for each galaxy are listed in Ta- 
ble ^ Col. (3) specifies the mass fractions (per cent) of the 
two components. Note that less than 10 per cent is required in 
all cases. The metallicities, ages, indices (H/3, Fe5335, Mgb), 
and colors {B-K, U-V) of the two populations are given in 
Cols. (4) -(10). The fuel Fuv (in Mq) in late evolutionary 
phases required to reproduce the observed UV-upturn is given 
in Col. (11). 

In Table || we present the relative contributions to H/3 and 
to the total fluxes at the wavelengths 5550 (V-band), 3600 (U- 
band), 2500, and 1500 A. The weight of the metal-poor com- 
ponent increases with decreasing A, with a maximum contribu- 
tion in the mid-UV (~ 2500 A) as the rising branch shortward 
2000 A is produced by the metal-rich component. It is impor- 
tant to emphasize that in the model proposed here, the Balmer 
line strengths and the UV upturn do not have a common origin: 
we obtain strong H/3 with metal-poor (hot) horizontal branch 
stars, while blue 1500-y colors are due to extremely hot old 
metal-rich stars. The trend of decreasing H/3 with increasing 
Mg2 and decreasing 1500-y (Burstein et al. 1988) is therefore 
not affected. 

In Table || we compare the resulting indices and colors of the 
composed models with the observed values. The indices H/3, 
Fe5335 and the colors B-K, U-V are in good agreement with 
observations. 

The fits to the observed SEDs are shown in Fig. ^ as thick 
solid lines, the dotted lines are the SEDs of the metal-poor 
and the metal-rich components. The spectra in the UV (open 
symbols) are from Burstein et al. (1988) and L. Buson (1999, 
private communication), the spectro-photometry in the optical 
(thin solid lines) is from D. Hamilton (1999, private communi- 
cation). In all cases, the observed spectra are well reproduced in 
the full wavelength range A = 1200-6000 A. Thus the amount 
of metal-poor stars required to reproduce the observed H/3 (see 
Table |l]) is consistent with the fluxes in the mid-UV, the model 
spectra perfectly match the minimum around 2500 A. In the 
following we briefly discuss the individual galaxies. 



3.2.1. M31 



As shown in Table |I] we reproduce the observables of M3 1 
(for the colors see NGC 4472) with Zsub = 0.005 Zq contribut- 
ing 3 per cent to the total mass of the population, which cor- 
responds to 6 per cent in V, 9 per cent at H/3 and 61 per cent 
at A = 2500 A. Wu et al. (1980), instead, claim that not more 
than 50 per cent contribution from a metal-poor component at 
A = 2500 A is compatible with the mid-UV flux of M31. The 
reason for this ostensible discrepancy is that the authors fix the 
metallicity of the metal-rich component to Z^aj = 0.5 Zq, which 
is lower than the value adopted here by more than a factor of 3 
(see Table |l]). As the SED of a more metal-rich population has 
lower flux in the mid-UV, our composite model requires a larger 
contribution from the metal-poor component. Using exactly the 
same prescriptions as these authors, we are consistent with their 
result. Due to the low metallicity of the metal-rich component, 
however, the model of Wu et al. (1980) fails in reproducing the 
other observational constraints of M3 1 : it results in colors that 
are too blue {B-K = 3.48, U-V w 1.21), Fe indices that are 
too low (Fe5335 « 2.01), and Balmer line strengths that are too 
high (H/3 = 2.07). 



3.2.2. NGC 4472 

The Virgo elliptical NGC 4472 has virtually the same prop- 
erties as M31, wherefore it can be modeled with the same set of 
populations. The model given in the second line of Table |l] is an 
alternative option for both objects assuming a higher contribu- 
tion from the metal-poor subcomponent and a higher metallic- 
ity of the major population. As a consequence, the contribution 
from the low-metallicity population increases to 70 per cent at 
2500 A and to 12 per cent in V. Wu et al. (1980) derive a lower 
contribution in V relative to 2500 A (50 per cent at 2500 A 
and 6 per cent in V), because the authors chose the globular 
cluster Ml 3 (NGC 6205) as a representative of the metal -poor 
population. As discussed in Section 2.2, this particular globu- 
lar cluster has an extreme blue horizontal branch morphology 
and exhibits very high fluxes in the mid-UV (see Fig. which 
implies a lower V over 2500 A flux ratio. 
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Table 1 
Model parameters 



galaxy 


comp. 


mass (%) 


Z(Zq) 


t (Gyr) 


H/3 


Fe5335 


Mgb 


B-K 


u-y 


Fuv (Mq) 


M31 


major 


97 


1.750 


12 


1.51 


3.20 


4.38 


4.34 


1.73 


0.003 




sub 


3 


0.005 


15 


3.09 


0.53 


0.78 


2.64 


0.53 




N4472 


major 


94 


2.000 


12 


1.48 


3.32 


4.54 


4.42 


1.81 


0.003 




sub 


6 


0.030 


15 


2.33 


1.02 


1.38 


2.91 


0.72 




N3379 


major 


99 


1.000 


14 


1.56 


2.82 


3.90 


4.06 


1.59 


0.002 




sub 


1 


0.005 


15 


3.09 


0.53 


0.78 


2.64 


0.53 




N4649 


major 


90 


1.550 


15 


1.38 


3.21 


4.41 


4.36 


1.73 


0.020 




sub 


10 


0.300 


15 


1.69 


2.08 


2.87 


3.53 


1.15 





Table 2 
Flux contributions 



galaxy 


H/3 


5550A 


3600A 


2500A 


1500A 


M31 


9 


6 


17 


61 


28 


N4472 


16 


12 


30 


70 


16 


N3379 


4 


2 


5 


33 


14 


N4649 


19 


17 


27 


16 






Note. — The Table gives the relative contributions from 
the subcomponent in per cent to H/3 (Col. (2)) and to the 
flux at the wavelengths specified in Cols. (3) -(6). 



Table 3 
Results 







H/3 




Fe5335 




B-K 


U- 


-y 


galaxy 


model 


observed 


model 


observedd 


model 


observed 


model 


observed 


M31 


1.65 


1.67 ±0.07 


3.04 


2.99 ±0.06 


4.26 




1.61 




N4472 


1.62 


1.62 ±0.06 


3.04 


2.84 ±0.08 


4.28 


4.28 


1.59 


1.57 


N3379 


1.62 


1.62 ±0.05 


2.77 


2.73 ±0.04 


4.05 


4.08 


1.55 


1.52 


N4649 


1.44 


1.40 ±0.05 


3.01 


3.01 ±0.05 


4.24 


4.26 


1.61 


1.61 



Note. — Indices and colors of the models given in Table |l]. Index data {la errors) are from Gonzalez 
(1993), B-K from Pahre (1998), and U-V from Trager (1998). 
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Fig. 5. — Spectral energy distributions as a function of wavelength. lUE data for A < 3300 A (circles) are from Burstein et al. 1988 
and L. Buson (1999, private communication); the spectro-photometry for A > 3440 A (thin solid lines) is from D. Hamilton (1999, 
private communication). The models from Table |l] are shown as thick solid lines, dotted lines are the SEDs of the metal-poor and the 
metal-rich components (100 per cent), respectively. 



3.2.3. NGC3379 

NGC 3379 exhibits equally high H/3, but bluer colors and 
lower metallic indices. The best fitting model thus has a lower 
•Zmaj = Zq and requires only 1 per cent of Zjub = 0.005 Zq. 

3.2.4. NGC 4649 

NGC 4649 has the sti'ongest UV upturn (bluest ISOO-V), re- 
quiring a factor of 10 more fuel at late evolutionary phases 
than the other cases (see column (1 1) in Table |l]). Given the low 
H/3 line strength, the best fit is obtained with a model contain- 
ing only old metal-rich populations (Zsub = 0.3 Zq). This object 
strongly supports the view that the UV rising branch and H/3 
line strength are not produced by the same kind of stars. 

3.3. Models for H(3>2 A 

The model grid in Fig. ^ extends to H/3 « 2 A, while some 
low-mass field ellipticals exhibit stronger Balmer lines. In prin- 
ciple such high H/3 can be obtained by increasing the weight 



of the low-metallicity component. For instance, the com- 
pact dwarf elliptical M32 (NGC 221) with H/3 = 2.37 ±0.12, 
Fe5335 = 2.54 ±0.20 (Kormendy & Bender 1999) can be mod- 
eled by the 14-Gyr old two-component population with Z^aj = 
2 - 3 Z0 and 20 per cent of Z^uh = 0.005 Zq. The SED of this 
model with such a large fraction of very metal-poor stars, how- 
ever is not compatible with the fluxes observed in the mid-UV 
(see also Burstein et al. 1984; Rose & Deng 1999). On the other 
hand, from a deep HST color magnitude diagram, Grillmair et 
al. (1996) conclude that there is no evidence for an intermedi- 
ate age population (f < 2 Gyr). It is worth mentioning, however, 
that M32 is a peculiar object that does not have the properties 
of the average elliptical galaxy population. 

Alternatively to the addition of metal-poor stars, metal-rich 
populations with blue horizontal branch morphologies could re- 
produce strong H/3 without invoking young ages. Such hot hor- 
izontal branch stars are indeed observed in the metal-rich glob- 
ular clusters NGC 6388 and NGC 6441 of the galactic bulge 
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(Rich et al. 1997) and in M32 (Brown et al. 2000). A possi- 
ble mechanism to produce these temperatures is enhanced mass 
loss along the RGB evolutionary phase. The investigation of 
this channel to obtain strong H/3 in old metal-rich populations 
will be the subject of a future paper. 

4. DISCUSSION 
4. 1 . Metallicity distributions in ellipticals 

The principle idea of this paper is to enhance the Balmer line 
strengths of old metal-rich stellar populations with a small frac- 
tion of metal-poor stars. The real metallicity distribution of the 
stellar populations in elliptical galaxies is difficult to assess ob- 
servationally, as the stars cannot be resolved. The HST color 
magnitude diagram of M32 (Grillmair et al. 1996) and the spec- 
troscopy of K giants in the Bulge (Rich 1988) show that these 
spheroids contain a tail of low metallicity stars. The closest gi- 
ant elliptical for which color magnitude diagrams are available 
is NGC 5128 (Harris, Harris & Poole 1999). Analyzing deep 
HST images in the outer halo of NGC 5128, the authors find 
a differential metallicity distribution that is well reproduced by 
two closed-box-like chemical enrichment scenarios implying a 
large number of extremely metal-poor stars. In particular the 
differential shape of the distribution implies a formation picture 
in which a metal-poor and a metal-rich population are formed 
separately from each other. 

This picture gets further support from a number of stud- 
ies that discover bimodal color and metallicity distributions of 
globular clusters in at least half of the early-type galaxy popu- 
lation (Zepf & Ashman 1993; Gebhardt & Kissler-Patig 1999). 
Spectroscopic and photometric investigations indicate that both 
populations are old (Cohen, Blakeslee & Ryzhov 1998; Kissler- 
Patig et al. 1998; Kissler-Patig, Forbes & Minniti 1998; Kundu 
et al. 1999; Puzia et al. 1999). More specifically, Puzia et al. 
(1999) find that the two populations in NGC 4472 are old and 
coeval, with metallicities ~ 0.05 Zq and ^ Zq. In the frame- 
work of hierarchical structure formation galaxies are built by 
mergers of smaller objects (e.g.. White & Rees 1978; White & 
Frenk 1991; Kauffmann, White & Guiderdoni 1993). A merger 
of coeval systems without newly induced star formation would 
result in a coeval composite population. The accretion of a 
dwarf galaxy by a larger system (minor merger) qualitatively 
explains the existence of a metal-poor subpopulation. 

There is an additional important effect that is independent 
of the assumed scheme of galaxy formation. As discussed by 
Greggio (1997), a sizable fraction of metal-poor stars from the 
outer parts of a galaxy contaminate the light coming from the 
center due to projection and orbital mixing (Ciotti, Stiavelli & 
Braccesi 1995). The authors find that projection effects lower 
the actual central metallicity by roughly 10 per cent. The metal- 
poor components in the composite models of M31, NGC 4472 
and NGC 3379 (Table |l]) dilute the metallicities of the metal- 
rich populations by 3, 6, and 4 per cent, respectively. Projec- 
tion effects alone thus may be sufficient to explain the amount 
of metal-poor stars considered in the present analysis. 

4.2. Evidence against recent star formation 

There are observational indications from the Fundamental 
Plane (Djorgovski & Davis 1987; Dressier et al. 1987; Bender, 
Burstein & Faber 1992, 1993; Renzini & Ciotti 1993) and the 
color-magnitude relation (Bower, Lucey & Ellis 1992) that limit 
the fraction of the young population to less than 10 per cent. 



As late star formation leads to low a/Fe ratios (Thomas, Greg- 
gio & Bender 1999; Thomas 1999), the addition of a young 
population provokes inconsistencies with the high a/Fe ratios 
observed in elliptical galaxies (e.g., Worthey et al. 1992). Also 
the redshift evolution of colors (Aragon Salamanca et al. 1993), 
of the Mg-cr relation (Bender, Ziegler & Bruzual 1996), of the 
Kormendy relation (Ziegler et al. 1999), and of the color gra- 
dients in cluster ellipticals (Saglia et al. 2000) leave little space 
for a significant contribution from a young subcomponent. The 
addition of an old, metal-poor population as discussed in this 
paper can reconcile these findings with strong Balmer lines. 

Disturbed field ellipticals that have strong H/3 absorption 
(Schweizer et al. 1990) and blue optical colors (Schweizer & 
Seitzer 1992) do not show any signature of recent star formation 
activity in the infrai-ed colors (Silva & Bothun 1998a, 1998b). 
As discussed by these authors, recent mergers may not have 
been accompanied by significant star formation, but metallicity 
effects are favored to explain the enhanced H/? line strengths. 

5. SUMMARY 

In this paper we show that strong Balmer lines can be pro- 
duced by composite populations that contain a small fraction of 
old metal-poor stars. 

The key to this result is that at low metallicities and old ages 
our SSP models show a steep increase of H/3 with decreasing 
metallicity, in accordance to what is observed in metal-poor 
globular clusters. We compute composite stellar population 
models that consist of an underlying old metal-rich population 
and a small fraction of an old metal-poor subcomponent. As- 
suming a 10 per cent contribution from metal -poor stars, we 
construct a grid of models in the H/9-Fe5335 and the \{j3-B-K 
planes that reproduce H/3 up to 2 A and covers most of the 
elliptical galaxy data. It should be emphasized that the same 
conclusion holds for all Balmer lines, e.g. H7 and H(5. More 
specifically, the data of all cluster and luminous field ellipticals 
can be explained by our models without invoking young ages. 
The scatter in Balmer line strengths is then caused by a spread 
in the metallicity and/or the weight of the metal-poor compo- 
nent. Most ellipticals can be modeled with less than 10 per cent 
of metal-poor stars as the data scatter about H/3 '-^ 1 .6 A. 

We further show that the spectral energy distributions of the 
composite models that reproduce the indices and colors of rep- 
resentative examples agree well with the observed spectra in 
the wavelength range 1200-6000 A. Our models are perfectly 
compatible with the observed minimum at 2500 A, that tightly 
constrains the possible amount of metal-poor stars. 

We conclude that the age-metallicity degeneracy for complex 
systems like elliptical galaxies still remains to be solved. The 
perturbation with old metal-poor stars to obtain strong Balmer 
lines is alternative to the addition of a young metal-rich popula- 
tion (e.g., de Jong & Davies 1997; Kuntschner 2000; Longhetti 
et al. 2000). The key to discriminate between the two options 
lies in the evolution of Balmer lines with redshift. A young 
population of 2 - 5 Gyr should leave its fingerprints with a sig- 
nificant peak of Balmer line strengths at redshifts z ^ 0.15-0.4, 
depending on the cosmology, which has not been detected so far 
(Ziegler & Bender 1997). The model presented here, instead, 
predicts Balmer lines to become monotonically stronger with 
redshift. 

We are grateful to S. Cassisi, D. Hamilton, D. Mehlert, and 
L. Buson, for providing electronic versions of models and data. 
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